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ABSTRACT 

Planetesimal formation is a common by-product of the star formation pro- 
cess. Taking the dynamical history of the solar system as a guideline - in which 
the planetesimal belts were heavily depleted due to gravitational perturbation 
with the giant planets - and assuming similar processes have take place in other 
planetary systems, one would expect the interstellar space to be filled with extra- 
solar planetesimals. However, not a single one of these objects has been detected 
so far entering the solar system, even though it would clearly be distinguishable 
from a solar system comet due to its highly hyperbolic orbit. The Large Synoptic 
Survey Telescope will provide wide coverage maps of the sky to a very high sen- 
sitivity, ideal to detect moving objects like comets, both active and inactive. In 
anticipation of these observations, we estimate how many inactive " interstellar 
comets" might be detected during the duration of the survey. The calculation 
takes into account estimates (from observations and models) of the number den- 
sity of stars, the amount of solids available to form planetesimals, the frequency 
of planet and planetesimal formation, the efficiency of planetesimal ejection, and 
the possible size distribution of these small bodies. 
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1. INTRODUCTION 

The study of protoplanetary and debris disks indicates that planetesimal formation is 
a common by-product of the star formation process, and there is observational evidence 
that in some cases it has led to the formation of planets. Taking the dynamical history 
of the solar system as a guideline, in which efficient dynamical ejection events triggered by 
gravitational perturbation with the giant planets were likely responsible for the depletion of 
the asteroid belt and the Kuiper belt (KB), and the formation of the Gort Cloud (OC), and 
assuming similar processes have taken place in other planetary systems, one would expect 
the interstellar space to be filled with extra-solar planetesimals. However, not a single one 
of these extra-solar planetesimals has been detected so far entering the solar system, even 
though it would be clearly distinguishable from a solar system comet due to its hyperbolic 
orbit resulting from the velocity of the Sun with respect to the Local Standard of Rest 
(LSR; vlsr — 16.5 km s~^). McGlynn & Chapman (1989) addressed this lack of "extra-solar 
comet" detections by comparing it to the expected number of detections based on a simple 
model that considers that stars have a spatial density of 0.1 pc~^, each of them ejects ~ 10^^ 
planetesimals into interstellar space, and the entire ensemble of planetesimals has a velocity 
distribution with respect to the Sun similar to that of the stars in the local neighborhood. 
Including corrections for the observability of these extra-solar comets, McGlynn & Chapman 
(1989) concluded that there is a missing extra-solar comet problem because a significant 
number of extra-solar comets should have already been detected. However, Jewitt (2003) 
argued that if each star ejects 10^^ planetesimals with r > 1 km (resulting in a number 
density of planetesimals of 10~^ AU~^) and a velocity relative to the Sun of 20 km s~^, one 
would expect 0.3 planetesimals per year to cross within 5 AU of the Sun (5 AU being the 
detection distance with Pan-STARRS for a low albedo comet-like body). Because a 1 km size 
object at 5 AU corresponds to about 24 mag. Jewitt (2003) concluded that it is no surprise 
none have been found so far and that Pan-STARSS would either detect a few interstellar 
comets over its 10 year lifetime or place limits on their number density. 

The last few years have witnessed a revolution in the fields planetary system studies 
(solar and extra-solar) with the completion of extended surveys of protoplanetary and debris 
disks, and most importantly, with the discovery and characterization of extra-solar planets 
and of Kuiper Belt objects (KBGs), shedding new light into the amount of solid material that 
might be available to form planetesimals, the frequency of planetesimal and planet formation, 
and the size distribution of small solar system bodies, respectively. These are the key features 
that affect the model calculations in McGlynn & Chapman (1989) and Jewitt (2003) in a 
fundamental way. In the future, cometary observations will witness a revolution on its own 
when Pan-STARRS and the Large Synoptic Survey Telescope (LSST) become operative, as 
they will provide wide coverage maps of the sky to a very high sensitivity ideal to detect 
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moving objects like comets (active and inactive). Because of these new developments, there 
is the need to revisit the interstellar comet estimates. 

To calculate the number density of planetesimals in interstellar space, McGlynn & Chap- 
man (1989) and Jewitt (2003) assumed a number density of stars that is independent of stel- 
lar mass, a fixed number of planetesimals each star contributes (also independent of stellar 
mass), and that all planetesimals are km-sized. We improve these approximations by taking 
into account the following observational and modeling results. 

• The number density of stars depends on the stellar mass (Section 12. ip . 

• Protoplanetary disks studies indicate that the amount of solids available to form plan- 
etesimals also depends on the stellar mass (Section l2.2p . 

• Debris disks surveys confirm that for a fraction of stars this solid material has accreted 
to form planetesimals. But debris disk surveys are limited by sensitivity to ~100 times 
the amount of dust in the solar system; this observational limitation, together with 
theoretical expectations, indicate that it is reasonable to assume all systems harboring 
giant planets also formed planetesimals (Section 12. 3p . 

• Giant planets eject planetesimals into interstellar space very efficiently (Section 12.41) . 

• Radial velocity surveys indicate that a fraction of stars harbor giant planets. Because 
the presence of massive planets seems necessary to provide a mechanism to scatter 
planetesimals into interstellar space and, as we mentioned above, it is reasonable to 
assume all systems harboring giant planets also formed planetesimals, we will assume 
that the stars contributing to the population of planetesimals in interstellar space are 
those that harbor giant planets (Section 12. 5p . 

• Not all planetesimals ejected are km-sized. We estimate the number density of planetes- 
imals with r > R based on the observed size distribution of the small body population 
in the solar system and on theoretical considerations (Section [3]). 

• Assuming the planetesimals are distributed isotropically in interstellar space and that 
they have an albedo of 6% (similar to that of inactive comets), we calculate the lu- 
minosity function of extra-solar planetesimals as observed from the Earth (Section 
i). 

• Finally, using the luminosity function above and simple considerations about the char- 
acteristics of the LSST survey, we estimate the expected number of "interstellar comet" 
detections (Section [5]). 
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2. MASS DENSITY OF EXTRA-SOLAR MATERIAL 

2.1. Number Density of Stars 

We adopt the stellar density calculated by Kroupa, Tout & Gilmore (1993): the number 
density of stars per pc^, out to ~ 130 pc from the Sun, within the midplane of the galaxy 
and with stellar masses between M* and M*-/-o?M* (in units of Mq), is given by 

n(M,) = i{M^)dM^ with, 
C(M,) = 0.035M-^-^ if 0.08 < M, < 0.5 
e(M,) = 0.019M-2-2 if 0.5 < M, < 1.0 
i{M^) = 0.019M-2-^ if 1.0 < M, < 100. 

The binary fraction of stars is approximately 50%. 

2.2. Total Mass Available to form Solids per Star 

The contribution of each star to the number of planetesimals in interstellar space is going 
to depend on the amount of solid material (primordial dust) that surrounds each star and 
that might be available to accrete into larger bodies. Andrews & Williams (2007) conduced 
a submillimeter survey of 170 Class II sources (pre-main-sequence stars with evidence of 
thick protoplanetary disks) in two star-forming regions, Taurus (<age>~l Myr) and p-Oph 
(<age>~0.7 Myr). It was found that the disk mass distributions are very similar with 
<Mdisk >~ 0.005 Mq and <Mdisk/M^ >~ 0.01 (where M^isk includes both gas and dust), 
for stellar masses ranging from 0.1 to 3 Mq (peaking at 0.4 Mq), and assuming a disk size 
of 100 AU and a gas-to-dust ratio of 100:1 (Andrews & Williams 2007). This is comparable 
to the minimum-mass solar nebula (~0.015 Mq - including gas and dust), which is the total 
disk mass required to account for the condensed material in the solar system planets. Even 
though Mdisk/M* decreases with stellar age from a few percent at 0.1 Myr to 30 x smaller 
value at 10 Myr, these Class II disk masses indicate that ~ 10~^M* of sohds could have 
been available to form planetesimals. Andrews & Williams (2007) note that their disk mass 
estimates adopted an opacity characteristic of mm-sized bodies. If larger bodies were to be 
present (cm or meter-sized), the opacity would be smaller than assumed and therefore the 
disk mass would be underestimated by up to an order of magnitude. Because the presence 
of these larger grains is currently unconstrained by the observations (except in the case of 
TW Hydra - Wilner et al. 2005), in this paper we will assume MgoUds ^ 10~^M*. 
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2.3. Fraction of Stars that Formed Planetesimals 

The disks of dust observeci^ around mature main-sequence stars of a wide range of 
spectral types indicate that, in some cases, there is evidence that the primordial dust particles 
in the protoplanetary disks have accreted to form planetesimals. This is inferred from the fact 
that the lifetime of the dust particles in these dust disk - under the effect of mutual collisions, 
Poynting-Robertson drag and radiation pressure - is 10^-10^ years, much shorter than the 
age of the star, 10''-10^° years, implying that the dust observed in these mature systems 
is not primordial, but must be replenished by a reservoir of dust-producing planetesimals. 
This is why these dust disks are known as debris disks. The characteristic dust temperatures 
inferred from the spectral energy distributions, together with the images of the two dozen 
of spatially resolved debris disks known to date, indicate that the debris disks have sizes 
comparable to the distribution of the small bodies in the solar system (from tens to few 
hundreds of AU). 

The Spitzer Space Telescope carried out extensive surveys to characterize the frequency 
and properties of the debris disks around more than 700 stars of different spectral types, ages 
and environment. For the purposes of this paper, we are interested in the fraction of stars 
that show evidence of dust excess emission, i.e., the fraction of stars that show evidence of 
harboring (dust-producing) planetesimal. 

• solar-type stars: for stars younger than 300 Myr, the disk frequency at 24 /im - tracing 
dust in the asteroid belt-like region - is 14.7%, decreasing to 2.7% for older stars; at 
70 /im - tracing dust in the KB-like region - the disk frequency is 6%-10% and does 
not change significantly with age (from a survey of 328 F5-K7 stars - Hillenbrand et 
al. 2008; Meyer et al. 2008; Carpenter et al. 2009). A different survey of 225 F0-G9 
stars older than 600 Myr yield disk frequencies of 4.2l^ 2% 24 /jm and 16.4l|®% at 
70 /im (Trilling et al. 2008). 

• A-type stars: for stars younger than 30 Myr, the disk frequencies are 42% at 24 fim - 
tracing dust at <5-50 AU - and 48% at 70 /im - tracing dust at 50-200 AU; for stars 
older than 400 Myr, these numbers decrease to 2% at 24 fim and 12% at 70 /xm (from 
a survey of 160 B9-A7 single stars - Su et al. 2006; Rieke et al. 2005). 

• Low-mass stars: the fraction of stars with debris disks decreases steeply for stars later 
than K2, with an excess rate for old M stars of 0% with upper limits (binomial errors) 



^Most of the 300 debris disks known to date are spatially unresolved and were identified from the presence 
of an infrared excess in the spectral energy distribution, attributed to the thermal emission from dust particles 
heated by the central star. 
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of 2.9% at 24 /im and 12% at 70 /iin (Gautier et al. 2007; Beichman et al. 2006; 
Trilling et al. 2008). This may be an observational bias because debris disks around 
this low-mass stars would have a peak emission at A > 70 fim and therefore remain 
too cold to be detected by Spitzer. 

• Binary stars: the frequency of debris disks around binary stars is similar or even higher 
than around single stars of similar type (from a survey of 69 A3-F8 binary systems - 
Trilling et al. 2007). 

The above debris disk fractions are lower limits because of the limited sensitivity of 
the Spitzer observations: Bryden et al. (2006) found that the frequency of dust detection 
increases steeply as smaller fractional luminosities (Ldust/L^) are considered, going from 
nearly 0% for Ldust/U = 10-^ to 2±2% for Ldust/U = 10"^ and 13±5% for Ldust/U = 
10~^; using the cumulative distribution of fractional luminosities and assuming that the 
distribution is a Gaussian, they concluded that the luminosity of the solar system dust (with 
Ldust/L^: = 10^'^ to 10~^) is consistent with being 10 x brighter or fainter than an average 
solar- type star, i.e., even though no solar system debris disks analog have been found due to 
limited sensitivity of the detectors, the observations indicate that debris disks at the solar 
system level might be common. 

Other observational evidence also indicate that planetesimal formation is a robust pro- 
cess that can take place under a wide range of conditions: debris disks are present around 
stars with more than 2 orders of magnitude difference in stellar luminosity, also in sys- 
tems with and without binary companions (as circumstellar or circumbinary disks), and 
around stars with a wide range of metallicities (unlike planet bearing stars that are strongly 
correlated with high stellar metallicities). This indicates that planetary systems harboring 
dust-producing planetesimals are more common than those with giant planet^, which would 
be in agreement with the core accretion models of planet formation where the planetesimals 
are the building blocks of planets and the conditions required for to form planetesimals are 
less restricted than those to form gas giants. 

Because of the above considerations and because, as we will see below, the mechanism 
of planetesimal ejection into interstellar space likely requires the presence of massive plan- 
ets, we will assume that the fraction of stars contributing to the population of interstellar 
planetesimals is determined by the fraction of stars harboring massive planets, rather than 
the fraction of stars with debris disks. 



^This is also in agreement with the lack of correlation between the presence of debris disks and the 
presence of planets (Moro-Martm et al. 2007) 
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2.4. Ejection of Planetesimals into Interstellar Space 

The dynamical history of the solar system can be relatively well constrained through 
models and observations and its highlights, with regard to the ejection of planetesimals, can 
be summarized as follows. 

• Before 10 Myr after the Sun was formed, while the Sun was still embedded in its mater- 
nal stellar cluster and before the gas in the primordial protoplanetary disk dispersed, 
Jupiter and Saturn formed and scattered the planetesimals in the Jupiter-Saturn re- 
gion out to large distances. Some of these planetesimals were completely ejected from 
the system after multiple close encounters with the planets. Others remained bound 
because their pcriccntcrs were lifted beyond the orbit of Saturn (avoiding subsequent 
encounters with the planets) due to gravitational perturbations from the stars and 
the gas in the stellar cluster. This population of bound planetesimals formed the OC. 
A model that accounts for the characteristics of the OC estimates that 75-85% of 
planetesimals in the Jupiter-Saturn region were completely ejected from the system 
(Brasser, Duncan & Levison 2006). 

• Between 10 and 100 Myr, O'Brien et al. (2007) estimate that the gravitational per- 
turbations from Jupiter and Saturn and the mutual perturbations amongst the largest 
asteroids depleted the asteroid belt by a factor of ~100, leaving behind an asteroid 
belt about 10-20 x more massive than today. 

• At ~ 700 Myr, the slow migration of the giant planets (due to the interaction of the 
planets with a massive trans- Neptunian planetesimal disk) forced Jupiter and Saturn 
to cross a mean motion resonance making their orbits unstable and triggering a quick 
rearrangement of the planets orbits into their current configuration (Gomes et al. 2005; 
Morbidelli et al. 2005; Tsiganis et al. 2005). During this process, secular resonances 
swept through the asteroid belt making their orbits unstable, scattering some asteroids 
into the inner solar system (producing the Late Heavy Bombardment - Strom et al. 
2005), while ejecting others out of the system, resulting in a depletion factor of ~10-20 
(i.e. depleting 90-95% of the asteroids that had remained). This planet rearrangement 
also resulted in the sudden outward migration of Neptune and its exterior mean motion 
resonances, that swept through the KB that as a consequence got heavily depleted. 

Overall, for the solar system it is estimated that only a very small fraction of negligible 
mass of the initial planetesimal disk was left behind. 

As can be inferred from the above description, the efficiency of planetesimal ejection 
is very sensitive to the dynamical history of the system, which in turn depends strongly on 
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the planetary configuration (the number of planets, their semimajor axis, masses and mass 
ratios) and the planetesimal disk (mass, extent and radial profile). The large diversity of 
planetary systems inferred from radial velocity observations warn us that there is no reason 
to assume that planetary systems share similar dynamical histories. However, the following 
considerations hint that the efficiency of planetesimal ejection may also be high for many 
other planetary systems. 

• The high occurrence of Jupiter-sized planets at small semimajor axis indicates that 
these planets have migrated significantly inward from their formation sites, a migration 
that likely resulted in the ejection of most of the planetesimals that interacted with 
the planet. 

• Highly eccentric planets are also common, however, there is no reason to expect plane- 
tary formation processes to yield high eccentricities and therefore it has been suggested 
that they result from periods of dynamical instability in multi-planet systems that lead 
to a sudden rearrangement of the orbits (Ford & Rasio 2008; Juric & Tremaine 2008). 
Models and observations of the solar system indicate that the sudden rearrangement 
of the giant planets orbits triggered the ejection a large number of planetesimals in the 
KB and the asteroid belt-regions. 

• Hints of transient events that may be due to dynamical instabilities are found in the 

observation of several anomalously massive debris disks around old stars with inferred 
levels of dust production that could not be sustained for the age of the system (Wyatt 
et al. 2007). 

For the purposes of this paper, we will assume that as in the solar system, most of the 
planetesimals in the planetary systems are ejected and only a small fraction of negligible 
mass is left behind. 



2.5. Fraction of Steirs with Massive Planets 

Because the mechanism of planetesimal ejection into interstellar space likely requires 
the presence of massive planets, and as we discussed above, we can assume planetesimals 
formed in all systems with planets, the fraction of stars contributing to the population of 
interstellar planetesimals is determined by the fraction of stars harboring massive planets. 

• solar-type stars: radial velocity searches indicate that 6.7% of solar-type stars have 
planets with masses > 0.1 Mj and semimajor axis < 5 AU (Marcy et al. 2005); from 
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the extrapolation of these surveys it is expected that ~19% harbor giant planets are 



• A-type stars: there are no firm statistics because radial velocity studies are complicated 
by the rotational broadening of the absorption lines, the decreased number of spectral 
features due to high surface temperature, and a large excess velocity resulting from 
inhomogeneities and pulsation. A way to circumvent these problems is by studying 
the cooler intermediate mass subgiants into which A-type stars evolve. Johnson et al. 
(2007a) found that for 1.3-1.9 M© subgiants (corresponding to F and A main-sequence 
stars) the frequency of planets with masses >0.8 Mj within 2.5 AU is 8.9±2.9% (= 
9/101), compared to 4.2±0.7% (= 34/803) for solar mass stars and 1.8±1% (= 3/169) 
for low mass-K and M dwarfs. These results indicate that A-F stars are five times 
more likely than M dwarfs to harbor a giant planet, although there is a possible bias 
from metallicity because the solar-type stars and the subgiants in the planet surveys 
are slightly metal-rich compared to stars in a volume-limited sample. After correcting 
for this effect, the above ratio decreases from 5 to 2.5 but is still significant (Johnson 
et al. 2007a). 

This positive correlation between the stellar mass and the occurrence of planets is 
in agreement with core accretion models for planet formation that predict that more 
massive stars result in more massive protoplanetary disks with the increased surface 
density of solid material in the midplane favoring the formation of planetesimals (Ida & 
Lin, 2005). Kennedy and Kenyon (2008) predict that the probability that a given star 
of mass 0.4-3 Mq harbors at least one giant planet increases linearly with stellar mass. 
For stars >3 the snow line quickly moves out to 10-15 AU before protoplanets 
form, limiting the formation of gas giant planets. 

Here we will assume that the fraction of A-type stars harboring planets is at least as 
high as for solar- type stars, ~ 20%. The effect of considering a higher planet fraction 
can easily be estimated from Table 1. 

• Low-mass stars: because ~ 77%±10% of stars are M-dwarfs, it is important to consider 
their possible contribution to the population of interstellar planetesimals. Johnson 
et al. (2007b) found that out of 300 M dwarfs monitored, only two have Jupiter- 
mass planets, while six have Neptune-Uranus mass planets (giving a total frequency 



^This comes from adopting the following probability for a star to have a planet of mass M at orbital 
period P: dN = CM"P^ dlnMdlnP, with a = -0.31±0.2 and /3 = 0.26±0.1. Assuming solar-type stars and 
planet masses O.SMjup < M < 15Mj„p, the frequencies are ~ 8.5% for a < 3 AU and ~19% for a < 20 AU 
(Gumming et al. 2008); if we were to consider larger semimajor axes, the frequencies would be ~ 23% for 
a < 40 AU and -33% for a < 100 AU. 
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of planets around M-dwarfs of ~ 8/300 ^ 2.7%). Even though the planet searches 
around M-dwarfs are limited by small number statistics, low-mass planets (Neptunes 
and super-Earths) appear to be relatively more common than Jupiter-mass planets 
(Bonfils et al. 2007; Endl et al. 2008). 

These observations agree with core accretion models that predict Jupiter-mass planets 
are relatively rare around M-dwarfs (due to the longer orbital periods, by the time 
the core reaches the critical mass of 10 Mq, the gas disk has dissipated and can no 
longer feed the planet's gaseous envelope - Laughlin et al. 2004; Ida and Lin 2005; 
Kennedy & Kenyon 2008), while Neptune-mass planets form more easily (Wetherill 
1996, Laughhn et al. 2004; Ida & Lin 2005 - the latter predicting a higher frequency 
of short-period ice giants around M-dwarfs compared to solar-type stars). 

For the purposes of this paper, we will assume that the fraction of M dwarfs with 0.1 
Mq< M < 0.6 Mq that can contribute to the population of interstellar planetesimals is 
~ 3%, and that because of the lower gravitational potential of M dwarfs, Neptune-mass 
planets can efficiently eject planetesimals. 

• Binary stars: the dependency of the frequency of planets on the presence of stellar 
companions is still under investigation. Bonavita and Desidera (2007) carried out a 
literature search of companions to the stars in Fischer and Valenti (2005) sample, 
with a uniform planet detectability for periods < 4 years. They found that: (1) the 
global frequency of planets in the binary sample (15/202 = 7.4%) is not statistically 
different from that of planets in the single star sample (5.3%). Even when making the 
conservative assumption that all stars without identified companions in the literature 
arc single, the frequency of planets in binaries is no more than a factor of 3 lower 
than that in single stars. (2) The frequency does not depend on the binary separation, 
except for close binaries. (3) For close binaries, the frequency of planets is lower than 
for single stars and wide binaries (however, the RV surveys are biased against binaries 
closer than ~2"). On going surveys, consisting on a radial velocity search for planets 
around spectroscopic binaries and an adaptive optics search for stellar companions to 
stars with and without planets, will soon shed light on the dependency of planets on 
binarity. 

In this paper, we will assume that only single stars contribute to the population of 
interstellar comets, adopting a binary fraction of 50%. 
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2.6. Expected Mass Density of Extra-Solar Planetesimals 

With tlie considerations discussed above, we can now calculate the contribution of stars 
of masses between and M^^+dM^ to the total mass density of extra-solar planetesimals. 
We assume that (1) the number density of stars per pc^ is given by Equation ([1]). (2) Each 
star harbored a protoplanetary disk with a total mass in solids of 10~^Mi, (Section I2.2p and 
all this solid material accreted into larger bodies (Section 12.31) . (3) 50% of the stars are 
binary, with 20% of the single A-K2 type stars and 3% of the single M-dwarfs harboring 
planets (Section [2.5p . able to eject most of the planetesimals out of the system (Section [2.4p . 

The results are listed in Table 1, adding up to a total mass density of interstellar material 
oimtotai = 2.2-10-^ Mo/pc^ = 4.5-1026 g/pd 



3. EXPECTED NUMBER DENSITY OF EXTRA-SOLAR 
PLANETESIMALS WITH r>R 

From the total mass density of extra-solar material {rritotai) derived above, we can esti- 
mate the number of planetesimals formed within a particular size bin by assuming a plan- 
etesimal bulk density of 1.5 g/cm^ (between that of pure ice and the density inferred for 
Pluto) and a size distribution. For the latter, we use the range of size distributions found for 
the small body population in the solar system as a guide of the ranges that may be found 
in the planetesimals ejected by other systems: 

• Asteroids. 

In the asteroid belt, the current size distribution is well known down to a size of 1 km, 
showing a wavy structure with enhancements at 3-4 km and 100 km, and a power-law 
index for the differential size distribution ranging from 2.25 to 3.8. Because this size 
distribution has evolved with time, more relevant is the size distribution at the time 
when the asteroids were ejected out of the system, i.e., at the time the solar system 
likely contributed to the population of extra-solar planetesimals. Using a collisional 
evolution model together with observational constraints, Bottke et al. (2005) concluded 
that the size distribution of this "primordial" asteroid belt followed a broken power 
law: n{r) oc r~'^^ ii r < rf, and n(r) oc r~'^^ if r > r^, where r is the planetesimal radius, 
Tb 50 km, qi 1.2 and q2 ~ 4.5. This distribution was established early on as a 
result of a period of collisional activity before Jupiter formed (few Myr), and a period 
of collisional activity triggered by the planetary embryos (10-100 Myr). 
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• Kuiper Belt Objects. 

The differential size distribution expected for KBOs, from coagulation models that 
take into account the coUisional evolution due to self-stirring, follows a broken power 
law with n(r) oc r^'^^ if r < ri, n(r) cc constant if ri < r < Tq, and n(r) oc r~'^^ if r 

> To, where r is the planetesimal radius and qi ^ 3.5 (resulting from the collisional 
cascade); for a fragmentation parameter, Qb > 10^ erg/g, q2 ~ 2.7-3.3, and ro ~ ri 
«i 1 km; while for Q^, < 10^ erg/g, q2 ~ 3.5-4, ri 0.1 km, and ro ~ 10-20 km (see 
review in Kenyon et al. 2008). Recent KBO surveys reviewed by Kenyon et al. (2008) 
yield q2 ^ 3.5-4 and r^ax ~ 300-500 km for the cold classical KB (a = 42-48 AU, 
perihelion > 37 AU, « 4°); g'2 ~ 3 and r^ax ~ 1000 km for the hot classical KB (i 

> 10°); and g2 ~ 3 and Vmax ~ 1000 km for the resonant population; in all cases, the 
transitional radius is Tq ~ 20-40 km (for albedo ~ 0.04-0.07). From a Hubble Space 
Telescope survey, Bernstein ct al. (2004) found = 2.9 and q2 > 5.85 for the classical 
KB, and qi < 2.8 and q2 = 4.3 for the excited KB, with < 50 km in both cases. A 
pencil-beam search for KBOs with Subaru found Vb ~ 45±15(p/0.04)^'^'^ km (where 
p is the albedo - Fuentes, George & Holman 2009). Finally, the most recent Subaru 
survey by Fraser and Kavelaars (2009), sensitive to KBOs with r > 10 km, found qi 
— 1.9, 52 = 4.8, and rb < 25-47 km (assuming a 6% albedo). 

Kenyon et al. (2008) argued that the comparison between the observational and the 
modeling results indicates that self-stirring alone cannot account for the observed size 
distribution in the KB and that dynamical perturbations have played a major role, 
suggesting that the size distribution in the KB was frozen after a major event of 
dynamical ejection. Fraser and Kavelaars (2009) argued that the large implies 
increased collisional evolution, while the large q2 suggests there was an early end of the 
accretion stage. Even though it is tempting to adopt the size distribution of the KBOs 
as a proxy of extra-solar planetesimals, the above indicates that the size distribution 
depends significantly on the dynamical/coUisional evolution of the population at the 
time of the ejection. 

• Comets. 

The size distribution of elliptic comets follows a power law of index q f« 2.9 for r > 
1.6 km, and if including cometary near-earth objects (thought to be extinct elliptic 
comets) q ^ 2.6. However, eUiptic comets are hkely collisional fragments of KBOs 
and therefore their size distribution is heavily evolved due to collisions and perihelion 
passages. More representative of the population of planetesimals ejected by the early 
solar system might be the less processed nearly isotropic comets, in particular those 
in their first perihelion passage. These comets are thought to originate from the OC 
that was populated in the early solar system by planetesimals in the vicinity of the 
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giant planets. Their size distribution is not known because of the small number of 
detections, although there are indications that it is shallower than the elliptic comets 
(see review by Lamy et al. 2004). 



(2) 



Based on the size distributions discussed above, and because it is not known the degree of 
dynamical/collisional evolution of the extra-solar planetesimals at the time they are ejected 
from other planetary systems, we consider a size distribution following a broken power-law 
with these range of parameters: 

n(r) oc r~'^^ if r < 
n{r) oc r~'^^ if r > 
gi = 2.0, 2.5, 3.0, 3.5 
q2 = 3, 3.5, 4, 4.5, 5 
Tfc = 3 km, 30 km, 90 km 
Tmax ~ 1000 km and r^m 1 /im 

(For simplicity, we use r^in = 1 /im for all stellar types, although for A-type stars the 
blow-out radius is ~10 /im.) 

Using the size distributions above, the expected number density of extra-solar planetes- 
imals with radius r > R {< rb) is 



4-92 r-fc ^ "J ' 4-131 L"" ^ rj, / J ^g-j 



1 ^ I b -"^ i' b ~ -I ^ I' max ' b \' b 

(see the Appendix A for details). 

We now use Equation ([3]) to calculate the total number density of planetesimals per 
pc^ with radius r > i? for the range of size distributions discussed above (in Equation ([2])). 
The results are shown in Figure 1 and Table 2, indicating that the number of planetesimals 
per pc^ with r > 1 km is in the range O(10^)-O(10^°), differing very significantly from the 
value used by McGlynn & Chapman (1989) of 10^"^ pc^'^. Using the latter number density, 
and including corrections for the observability of extra-solar comets, McGlynn & Chapman 
(1989) concluded that there is a missing extra-solar comet problem because a significant 
number of extra-solar comets should have already been detected. However, our estimate 
of the number density of planetesimals is several orders of magnitude smaller so it is not 
surprising that extra-solar comets have not been detected so far. 
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=4.5x1 O^^g/pc^ p=1.5g/cm^ 




0.1 1.0 10.0 

R[km] 



Fig. 1. — Total number density of planetesimals per pc^ with radius r > R calculated from 
Equation (fT6!) for the range of size distributions in Equation ([2]). The colors and line types 
correspond to qi = 2.0 (light blue), 2.5 (red), 3.0 (blue), and 3.5 (green); g2 = 3 (sohd), 3.5 
(dashed), 4 (dotted), 4.5 (dash-dotted), 5 (long dashed). The different panels correspond 
to different break radius: Vh = 3 km (top panel), r;, = 30 km (middle panel), and = 90 
km {bottom panel). In all cases, the total mass density of extra-solar material is rritotai = 
2.2-10^^ Mq/pc^ = 4.5-10^^ g/pc'^ (derived in Section [2T6|) . the planetesimal bulk density is 
p = 1.5 g/cm^ rmax = 1000 km, and r^i„ = 1 /im. 
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EXPECTED LUMINOSITY OF EXTRA-SOLAR PLANETESIMALS 



Using the number density of planetesimals in interstellar space derived in Section [3l 
we can calculate the distribution of luminosities as observed from the Earth (see details 
of the calculation in the Appendix B). This will allow us to estimate how many of these 
objects might be detected, given the limiting magnitude of the survey under consideration. 
We assume the planetesimals are distributed isotropically and that their size-magnitude of 
brightness relation is 

m = K + 2.5logio{a-^D-^), (4) 

where D is the planetesimal diameter in km, a is the heliocentric distanc^ in AU, and K 
= 18.4 mag, which assumes that all objects have an albedo of 6%. This albedo is typical 
of inactive comets, therefore, this size-magnitude of brightness relation would only be valid 
for inactive planetesimals (i.e., with no comae). Following Fraser & Kavelaars (2009), the 
cumulative luminosity function, i.e., the number of objects per deg^ brighter than a given 
magnitude (m < nimax) is given by 



with, 



6ilO— ^""^ +6210^ 



+ &3 



-K(qi-l) 



(gi-l)(5-2g0 



-K)/10\5-2qi 



a 



5-2gi 
1 



(5) 



(6) 



A 



10 



-K(q2-1) 



(g2-l)(5-2g2) 



a 



-K)/10\5-2q2 



(7) 



AD 



l-'?2 



qi - q2 



3(g2-l)(gi-l) 



(D^^IO^" 



-K)/10\3 



A 



1.5 ■ 10 

3 • 1018 



13\ 3 



f— ) (2-10-^)^ 
V180/ ^ ^ 



9(47rp) ^mtotai {Dt 



,(?2-4 



_3 l( Dmax '\4-g9 1 1 _|_ 3 M _ ( Drrmi\4:-qi] 

4-92 Lv Db ^ J 4-gi ^ Dt ^ J 



(9) 



where the planetesimal diameter D is in km, its heliocentric distance a is in AU, its bulk 
density p is in g/cm^, and the total mass density of planetesimals rritotai is in g/pc^. Equations 
&, (IZD> &, and Q> are only valid for qi ^ 1, gi ^ 2.5, q2 7^ 1, and q2 7^ 2.5. 



'The distance to the Earth is approximated as the hehocentric distance. 
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Figure 2 shows the number of planetesimals brighter than a given magnitude [m < 
nT'max) per deg^ calculated from Equations ([5]), ([6]), ([7]), ([8]), and Qj, with the parameter 
ranges in Equation ([2]). These are the cumulative luminosity distributions that would cor- 
respond to the cumulative size distributions in Figure 1. As mentioned above, we assume 
the planetesimals have an albedo of 6%, i.e., they are inactive. To fulfill this constrain we 
consider planetesimals at heliocentric distances beyond 5 AU, i.e., oi = 5 AU and a2 — > oo. 
Table 2 shows the expected number of planetesimals per deg^ with m < 24.5 (the limiting 
magnitude of LSST for a single visit), for different values of the minimum and maximum 
planetesimal diameter {Dmin and Dmax)- 



LSST is expected to achieve a depth of 24.5 magnitudes for each single visit. Table 2 
shows that for the broad range of size distributions considered, the expected number of extra- 
solar planetesimals per deg^ brighter than this limiting magnitude is O(10^^°)-O(10^^). In 
the total area that will be covered by LSST, 2-10'' deg^, we would expect to find O(10~^)- 
O(10~^) extra-solar planetesimals. But we need to consider the following: 

• LSST will observe the sky 1000 timeJ^ in 10 years. Considering that the velocity of the 
extra-solar planetesimals with respect to the Sun is approximately the velocity of the 
Sun with respect to the LSR, vlsr ~16.5 km s~^, we expect the number of extra-solar 
planetesimals inside 5 AU to "refresh" approximately once per year. The numbers in 
Table 2 were calculated assuming albedos of 6%, corresponding to inactive comets, and 
to be consistent with this the objects are assumed to be at distances beyond 5 AU. 
This means that the "refreshing" of the population of extra-solar planetesimals within 
the 10 year lifetime of the LSST survey increases the estimated number of detections 
by less than a factor of 10. 



• The gravitational focusing factor at a heliocentric distance R is Gf = -\ — ^;T^j ? 
so this effect can increase the number of expected detections by a factor of ~ Gj. For 
R = 5 AU, ~ 2; because the objects are at distances beyond 5 AU, gravitational 
focusing could increase the number of detections by less than a factor of 10. 



5. 



DISCUSSION 




Because the extra-solar planetesimals are moving objects, we are limited by the single visit magnitude 
and in this case the number of visits does not increase the depth. 



-17- 




Fig. 2. — Number of planetesimals brighter than a given magnitude (m < rrimax) per deg^ 
for the size distributions in Figure 1 (same color code). We assume the planetesimals have 
an albedo of 6%. 
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The above two considerations lead to an increase in the expected number of detections of 
less than a factor of 100, i.e. we expect LSST to find at most O(10~^)-O(10~^) extra-solar 
planetesimalj"] during the duration of its lifetime. 

The effect of the refreshing of the extra-solar planetesimal population mentioned above 
can be calculated explicitly by estimating the expected incoming flux of extra-solar planetes- 

imals larger than a given radius R. This flux is / n{r)a{r)vdr, where v is the planetesimal 

R 

velocity with respect to the Sun {v ~ vlsr ~16.5 km s~^), n{r) is the planetesimal differen- 
tial number density (from Equations (fT3|) and (fT5l) ). and a is a "cross-section", cr(r) = vra^. 
In this lalter expression, a(r) is the maximum distance at which a planetesimal of radius r 
could be detected, a(r) = (2r)^''^10*^™""'"'~^-'/^°, where rrimax is the limiting detectable mag- 
nitude, r is the planetesimal radius in km and K = 18.4 (assuming an albedo of 6% - see 
Section m and Equation (jlj)). For rrimax = 24.5 mag (limiting magnitude of LSST for a single 
visit), and using Equations ffT^ and (HM . we get that the number of incoming detectable 
extra-solar planetesimals larger than radius R, per year, is 



J n{r)a{r)vdr = 4.5 • 10 

R 
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-gi+2 p-gi+2 \ / -92+2 .,-92+2 

A/ I J_b jj^ ) _|_ /I ' ' ^ 

11 „ I o „ I o ; 2 



-gi + 2 -qi + 2J '\-q2 + 2 -gs + 2 / 

(10) 

where A[ and A'^ are given by Equation f[T^ (with p = 1.5 g/pc^, rritotai = 4.5-10^^ g/pc^ 
- see §2.61) . and R, and r^ax are in cm. The results, shown in Figure [3l indicate that 
the likelihood of detection is larger for small objects close to the Sun, compared to larger 
bodies at a wider range of distances; the faintness (and hence smaller a) of small objects is 
compensated by their larger number density. 

We conclude that the fact that not a single extra-solar comet has been detected so far 
is not in contradiction with current knowledge (as pointed out by Jewitt 2003), and that the 
probability that LSST will detect a single inactive extra-solar comet during the duration of 
its lifetime is very small (~ 0.01%-1%): a non-detection of inactive comets is expected from 
our current knowledge of the distribution of stellar masses, the typical protoplanetary disk 



^^The upper limit could be 10 times larger if we were to assume that all the solids in the protoplanetary 
disks accrete into objects with a very narrow range of sizes, from 1 to 5 km with a break radius of = 3 km, 
but this very narrow size distribution is not supported by observations of small bodies in the solar system, 
nor by core accretion models of planet formation that predict the formation of larger bodies. 
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Fig. 3. — Expected number of incoming detectable extra-solar planetesimals larger than 
radius R, per year, for the size distributions in Figure 1 (same color code). We assume the 
planetesimals have an albedo of 6%. 
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mass (assuming opacities typical of mm-sized particleJ^. the frequency of planetesimals and 
planets, and the distribution of sizes of small bodies. Note, however, that the likelihood of 
detection is significantly larger for small objects close to the Sun. At these small heliocentric 
distances, the expected number of small detectable objects in Figure [3] is underestimated 
because of increased gravitational focusing, and more importantly, because these objects 
might have become active and more easily detectable. The study of how cometary activity 
affects the expected number of detections is left for future work. 



5.1. Caveats 

In our calculation we have assumed that extra-solar planetesimals have albedos of 6%, 
similar to inactive comets. To be consistent with this assumption we have considered extra- 
solar planetesimals located beyond 5 AU, i.e., before comets start outgassing. However, 
the properties of extra-solar planetesimals are not known. In the solar system, albedos can 
vary from 2% up to 70%. For example. Centaurs have no comae, either because they have 
already lost their volatiles, or because their surfaces consist of non- volatile complex organics 
and silicate mantles processed by the bombardment of energetic particles (Jewitt 2004). 
Extra-solar planetesimals have long interstellar travel times and is expected that they are 
exposed to cosmic rays for long periods of time. It is therefore reasonable to assume that 
their crusts are highly processed and, as a consequence, have low albedos before they become 
active - as we have assumed in our calculations. 

But we do not know at what heliocentric distance extra-solar planetesimals activate and 
how bright they might become because the composition of the ices could be very different 
from solar system comets; for example, CO, if dominant, could enhance the comet activity 
far from the star. New nearly isotropic comets (i.e., those in their first perihelion passage) 
are anomalously bright, probably because of their high volatile content (accreted from the 
interstellar medium), and/or the blow-off of their cosmic-ray-processed crust, and/or the 
amorphous-to-crystalline H2O ice phase transition (occurring at 5 AU during their first 
perihelion passage - Dones et al. 2004). Extra-solar comets might also become unusually 
bright when they become active and their highly processed crust is blown off. This makes 
the estimate of how many active extra-solar comets might be detected very challenging. 

Another issue to keep in mind is that extra-solar planetesimals may break easily. From 
modeling and observations of inactive nearly isotropic comets, it is found that 99% become 



^^If larger particles were found to be common in the protoplanetary disks, the disk masses would increase 
affecting the estimated number of extra-solar comets. 
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unobservable as they evolve inward from the OC, probably due to breaking into smaller 
pieces; as these comets become inactive, they only have a 1% chance of becoming dormant, 
the other 99% must have disrupted; this 99% comes from comparing the number of returning 
comets to those predicted based on the number of new nearly isotropic comets (Levison et 
al. 2002; Dones et al. 2004). For comparison, Jupiter family comets do not seem to break 
so easily and about 60% become dormant]^ compared to 1% (Levison et al. 2002). If extra- 
solar comets are similar to the unaltered nearly isotropic (OC) comets, it may be possible 
that a high fraction of them break up and become unobservable before we have the chance 
to detect them. 



APPENDIX A 



DERIVATION OF THE NUMBER DENSITY 



Using the size distributions in Equation ([2]), we calculate the expected number density of 
extra-solar planetesimals with radius r > R. We adopt a broken power law for the mass 
distribution]^! 

n(m) = Ai'm~°'^ for m < ruh 

_ (11) 
n(m) = for m > m^, 

where Ai and A2 are calculated from continuity, AimV^^ = A2m7"'^, and from normaliz- 



ing to the expected mass density over the full mass range, rritotai = J Aim "^mdm + 
J A2m~°'^'mdm, yielding 



Al 



— an+2 — Qo+2 -a-1 +2 —a-, +2 

b -02+2 -ai+2 /-.^N 
. _ rritotai ^ ' 

^2 -ay+2 -00+2 -a-, +2 -ai +2 ' 



+ m 



-a2+2 ' b -Qi+2 



""^^Possible reasons are because they contain less volatiles (due to being stored at higher temperatures), 
and/or they are more porous (because of more frequent coUisions) and therefore less susceptible to disruption 
from the build-up of volatiles, and/or they are subject to smaller temperature gradients (Levison et al. 2002). 

^"^We start the discussion in terms of the mass distribution (rather than the size distribution) because 
we use the total mass density of extra-solar material {rritotai calculated in Section 12. 6p to normalize the 
distribution. 



-22- 



Assuming all planetesimals have the same bulk density, p, the above mass distribution is 
equivalent to a size distribution of 

nir) = A\r~^^ for r < tk 

(13) 

n(r) = for r > r?,, 

where A'^, qi and q2 are derived from n{m)dm = n{r)dr and given by, 

= 3°2(47rp)-"2+iA2, 



with gi = 3ai — 2 and g2 = 3^2 — 2. Substituting the expressions in Equation (fT2i) into ([T 

^Tcprf and = ^ 

9(47rp)-imtoi,;r,^^-' 



and using = Ivrprf and = we get 



_\(rmax \4:-g2 _ 11 J 2_ri 

Q 

/in 



9(47rp) ^rritotair] 



(15) 



^2 



.[(1^)4-92 + 



^4-gil 



The total number density of planetesimals with radius r > R {< ri,) is given by A^r>K 
J A[r~'^^dr + J A2r~'^^dr, and from the expressions in (fT5|) we get, 

9(47rp)-imtotai 



.[(IH^)4-92_1] + _J^^[1_ 



i — qi i- — q2 



APPENDIX B 
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DERIVATION OF THE LUMINOSITY FUNCTION 

The number of planetesimals with radius r to r + dr (with r in cm) per pc^ is given by 
Equations (fT3|l and (fT5|) . with the possible range of parameters listed in ([2]). If we assume 
the planetesimals are distributed isotropically, the number of plantesimals with radius r to 
r + dr (in cm) and heliocentric distance a to a + da (in AU) per deg^ would be given by 

n{r)dr f^^^j^r-j An a? da ^ /igoN^ ' where ^g^^^Qi" is the number of pc in an AU and An (^)^ is 

the number of deg^ is a sphere. For the derivation of the luminosity function, we need to have 
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the size distribution in terms of the planetesimal diameter D in km. The transformation is 
given by n{r)dr = S{D)dD^^^^; from Equations f|T3|) and f|T5l) . we get 



with, 



A'' 



S{D) = AID-"' for D<Di, 
S{D) = A'^D-'^' for D>Db 



Dmax ^ 4- 



-52 



-(2 • iQ-^Y^ 



( Dmax ^ 4- 



4-52 l-^ -Db 



'<?2 



— — -(2-10-^)^^ 



(17) 



The number of objects per deg^, with diameters from Dmm to Dmax (in km) and heho- 
centric distances from ai to 02 (in AU), is given by 
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(19) 
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iP-l 



A'i 



Djnax 



D-'i'dD+ j D-'^'dD 



^2 Dfnax 

da + j j D-^^dD 



where (in AU) is the maximum distance at which a planetesimal of diameter (in km) 
can be detected, determined by the hmiting magnitude of the survey. 

We now follow Eraser and Kavelaars (2009) to calculate the luminosity function corre- 
sponding to the size distribution in Equation (fT9l) . As in their calculation, we assume that the 
size-magnitude of brightness relation is m = if + 2.5/o(y'io(a~'^-D~^), where D is the planetesi- 
mal diameter in km, a is the heliocentric distanc in AU, and K = 18.4 mag, which assumes 
that all objects have an albedo of 6%. Using the relation above, the maximum distance at 
which a planetesimal of diameter (in km) can be detected is = /}y^io('"""'^~^)/^°. 



^■^The distance to the Earth is approximated as the hehocentric distance. 
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FoUowing Fraser & Kavelaars (2009), the cumulative luminosity function, i.e. the number of 
objects per deg^ brighter than a given magnitude (m < nimax) is given by. 



with, 



61 10^™'"'^^ + 6210^"™^ + 63, 
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(24) 



D-min \ 4 — 1 



where the planetesimal diameter D is in km, its heliocentric distance a is in AU, its bulk 
density p is in g/cm^ and the total mass density of planetesimals rritotai is in g/pc^. Equations 
dZH), (122D, ([23]) and ([21]) are only valid for gi 7^ 1, gi 7^ 2.5, 92 7^ 1, and q2 ^ 2.5. 
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Table 1. Expected Mass Density of Extra-Solar Material 



Spiype 


Mass Range 


/ e(M,)10-4M,dil4 


Fraction of 


Mass Density 




(Mo) 


(Mo/pc=^) 


Gontributing Stars 


(Mo/pc=^) 


M 


0.1-0.5 


2.1-10-*^ 


0.03x0.5 


3.1-10-^ 


MK2 


0.5-0.8 


6.1-10-^ 


0.03x0.5 


g.i-io-*^ 


K2G 


0.8-1 


4.3-10-^ 


0.2x0.5 


4.3-10-^ 


GF 


1-1.8 


9.1-10-^ 


0.2x0.5 


9.1-10-^ 


A 


1.8-2.9^ 


5.1-10-^ 


0.2x0.5 


5.1-10-^ 


AFGKM 


0.1 2.9 






2.2-10-^ 



^We do not consider larger stellar masses because for M^, >3 Mq the snow line move 
quickly beyond 10-15 AU before protoplanets form, limiting the formation of gas giant 
planets. 
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Table 2. Expected N{r > 1 km) and A^(m < 24.5) 





r, 


■nin 


(km) 


N{r > 1 km) 
(pc ^) 




A^(m < 24.5) 
(deg-2) 


3 km 


i 


/im 


1 nnn 
iUUU 


o 1 n7 p; /I 1 n9 
Z-iU — 0.4-iU 


o 
A 


R 1 n~io o /I 1 n^8 
o-iU — Z.4-iU 




1 
i 


/im 


ouu 


/I T n7 /I 1 nS 





i"iU ^Z.4"iU 




i 


/im 


ou 


4.i-iU — O.O-iU 





1 n^9 9/11 n— 8 




1 


/im 


5 


5.3-10 -8-10^ 


3 


7-10 -3.4-10 * 




1 


km 


1000 


2-10^-8.8-109 


2 


6-10-i°-3.3-10-* 




1 


km 


50 


4.1-10M-109 


5 


3-10-9-3.4-10~8 




1 


km 


5 


5.5-10M.3-10^° 


3 


9-10-*-7.7-10-^ 


30 km 


1 


/im 


1000 


2.4-10M.8-109 


1 


l-10-i°-8.3-10-9 




1 


/im 


500 


4.8-10M.8-109 


2 


3-10-i°-8.4-10^9 




1 


/im 


50 


6.7-10^-2.2-109 


2 


3-10-9-1. 0-10^^ 




1 


km 


1000 


2.4-106-2-109 


1 


l-10-i°-9. 5-10-9 




1 


km 


50 


6.7-10^-2.6-109 


2 


3-10-9-1. 2-10-s 


90 km 


1 


/im 


1000 


8.2-10M-109 


7 


1-10-11-5.0-10-9 




1 


/im 


500 


1.7-10M-109 


1 


5-10-1° -5.1-10-9 




1 


km 


1000 


8.2-10M.1-109 


7 


1-10-11-5.4-10-9 



^ N{r > 1 km) is the expected number of planetesimals per 
pc^ with radius r > 1km. N{m < 24.5) is the expected number 
of planetesimals per deg^ with magnitude brighter than 24.5 (the 
limiting LSST magnitude for a single visit). These values are 
calculated using Equation (fT6ll for N{r > 1 km) and Equations 
(EQ]), (lai), (1221), ([23]), and (Elj) for N{m < 24.5), for the range of 
size distributions discussed in Section [21 n{r) oc r~'^^ for r < Vj, 
and n(r) oc r~'^^ for r > r^, with qi = 2.0-3.5, q2 = 3-5, and Vh 
= 3km, 30 km, and 90 km. The total mass density of extra-solar 
material is mtotai = 2.2-10"'^ Mq/pc^ = 4.5-10^^ g/pc^ (derived 
in Section [22]) , the planetesimal bulk density is p = 1.5 g cm~'^ 
and the albedo is 6% (corresponding to inactive comets). 



